A Dual-Species Bose-Einstein Condensate of 87 Rb and 133 Cs 
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We report the formation of a dual-species Bose-Einstein condensate of 87 Rb and 133 Cs in the 
same trapping potential. Our method exploits the efficient sympathetic cooling of 133 Cs via elastic 
collisions with 87 Rb, initially in a magnetic quadrupole trap and subsequently in a levitated optical 
trap. The two condensates each contain up to 2 x 10 4 atoms and exhibit a striking phase separation, 
revealing the mixture to be immiscible due to strong repulsive interspecies interactions. Sacrificing 
all the 87 Rb during the cooling, we create single species 133 Cs condensates of up to 6 x 10 4 atoms. 

PACS numbers: 03.75.Mn, 03.75Hh, 05.30Jp 



Quantum degenerate mixtures [lHl2| of atomic gases 
are currently the subject of intensive research, exhibiting 
rich physics inaccessible in single species experiments. 
Phase separation of the constituent gases is a particu- 
larly dramatic example [H, lil-[lo|. In an optical lattice, 
the myriad of quantum phases for two different atomic 
species [HI, Q3] goes far beyond the seminal observations 
of the supernuid to Mott-insulator transition [Hj and of- 
fers many intriguing applications. Amongst these appli- 
cations, the creation of degenerate samples of heteronu- 
clear ground state molecules is of paramount interest [l6| . 
Possessing long-range, anisotropic dipole-dipole interac- 
tions, such molecules offer a wealth of new research av- 
enue s [171, i ncluding quantum simulation and computa- 
tion [13, [l9j. At present the most promising route to 
these riches exploits a mixed species quantum gas in 
which the constituent atoms are associated into ground 
state molecules [lo|. Remarkably, using a combination 
of magneto-association on a Feshbach resonance [2l| fol- 
lowed by optical transfer to the rovibrational ground 
state, this association can be achieved with near unit 
efficiency, as demonstrated using KRb Q and Cs2 0. 

A quantum degenerate mixture of 87 Rb and 133 Cs 
would provide an ideal gateway into the realm of dipo- 
lar molecular quantum gases. Both species have been 
separately condensed and the intraspecies two-body in- 
teractions are now well understood [24], [25| . Numerous 
interspecies Feshbach resonances have been observed , 
permitting, in principle, the control of the interspecies in- 
teractions and the creation of RbCs Feshbach molecules. 
Crucially, the RbCs molecule in the rovibrational ground 
state is stable against atom exchange reactions of the 
form RbCs + RbCs Rb 2 + Cs 2 , unlike the KRb 
molecule [I?! • However, associated with the rich Fesh- 
bach spectrum in 133 Cs are high two- and three-body 
inelastic collision rates which complicate the evaporative 
cooling. Indeed, to date Bose-Einstein condensation of 
133 Cs has required sophisticated laser cooling techniques, 
and has been achieved by only two groups [H, [29| . 

In this Letter we demonstrate the production of a 
dual-species Bose-Einstein condensate of 87 Rb and 133 Cs, 
thereby realizing the perfect starting conditions for the 
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FIG. 1: (Color online) (a) Geometry of the levitated crossed 
dipole trap used in the experiment, (b-d) Typical trapping 
potential employed towards the end of the evaporation se- 
quence (the power of both beams is 120 mW, the magnetic 
field gradient is 35 G/cm and the bias field is 22.4 G). The 
contour plot in (b) shows the potential for 87 Rb. The cut 
along z in (d) highlights the tilting of the trap. 



creation of a quantum gas of dipolar RbCs molecules. 
Our method uses simple established techniques and ex- 
ploits the efficient sympathetic cooling of 133 Cs via elastic 
collisions with 87 Rb, initially in a magnetic quadrupole 
trap and subsequently in a levitated optical trap. The 
method can easily be adapted to produce large single 
species 133 Cs condensates, offering an alternative route 
to quantum degeneracy for this atom. 

At the heart of our experiment is a levitated crossed 
dipole trap [30| located at the centre of an ultrahigh vac- 
uum (UHV) glass cell (Fig.[T]). Two 1550 nm beams de- 
rived from a 30 W fibre laser (IPG ELR-30-LP-SF) are fo- 
cussed to waists of ~ 60 fim and intersect ~ 80 fim below 
the unbiased field zero of a magnetic quadrupole poten- 
tial. The power in each beam is independently controlled 
using acousto-optic modulators. The magnetic potential 
generated by the quadrupole and bias coils allows atoms 
in the trap to be levitated against gravity. Both species 
can be levitated simultaneously due to the fortuitous co- 
incidence of the magnetic moment to mass ratios for the 
87 Rb \F = l,m F = +1) and 133 Cs |3,+3) states. When 
loading the trap the power in each beam is set to 6 W. 
This creates a trap 90 /xK (125 /zK) deep for 87 Rb ( 133 Cs). 
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The difference in the trap depths arises from the different 
polarizabilities at 1550 nm and makes the trap potential 
suitable for sympathetic cooling of 133 Cs. For this load- 
ing potential, the axial and radial trap frequencies for 
87 Rb ( 133 Cs) are v z = 135 Hz (127 Hz) and v r = 680 Hz 
(640 Hz) respectively. 

To load the dipole trap, we first collect ultracold mix- 
tures of 87 Rb and 133 Cs in a two-species magneto-optical 
trap (MOT) in the UHV cell Up to 3 x 10 8 87 Rb 
atoms are collected. The number of 133 Cs atoms in 
the MOT is stabilized between 5 x 10 5 and 3 x 10 8 by 
monitoring the MOT fluorescence and actively control- 
ling the repump light level. This allows precise control 
of the species ratio in the mixture. Following a com- 
pressed MOT and molasses phase the atoms are optically 
pumped into the 87 Rb |1,-1) and 133 Cs |3, -3) states 
and the mixture is loaded into a magnetic quadrupole 
trap. The trap gradient is adiabatically increased from 40 
to 187 G/cm, prior to using forced RF evaporation to cool 
the 87 Rb gas. To probe the mixture, the trap is switched 
off and resonant absorption images are captured for both 
species in each iteration of the experiment using a frame 
transfer CCD camera [31]. The trap depth set by the RF 
frequency is three times greater for 133 Cs than for 87 Rb 
and therefore allows the selective removal of hot 87 Rb 
atoms. As the 87 Rb gas is cooled we observe efficient 
sympathetic cooling of 133 Cs via interspecies elastic colli- 
sions [HI. The efficiency, 7, ranges from 5.2(4) to 11.3(4) 
depending on the initial number of 133 Cs atoms (the ef- 
ficiency is defined as 7 = — [^[n^/n] •> wnere Do (^0) 
and D (TV) are the initial and final phase-space densities 
(numbers) respectively). The high efficiency of the sym- 
pathetic cooling implies a large interspecies elastic colli- 
sion rate. This is consistent with the latest theoretical re- 
finement of the interspecies collision potential [HI which 
shows excellent agreement with the published Feshbach 



spectrum [26|, |34| and indicates a background interspecies 
scattering length of ~ 650 ao (ao = 0.0529 nm) for the 
states used in this work. We continue this cooling until 
Majorana losses begin to limit the lifetime of the sample. 
At this stage the crossed dipole trap is loaded via elastic 
collisions by simply reducing the magnetic field gradient 
to 29 G/cm, just below the 30.5 G/cm (31.1 G/cm) re- 
quired to levitate 87 Rb ( 133 Cs) [35]. At the same time the 
final RF frequency is adjusted to control the composition 
of the mixture loaded into the dipole trap. 

To create 133 Cs condensates a mixture containing 
2.1(1) x 10 8 87 Rb and 2.0(2) x 10 7 133 Cs atoms is loaded 
into the quadrupole trap. The evaporation efficiency for 
87 Rb alone in the magnetic trap is 2.6(2) (O in Fig. [2] 
(i)). For the two species case the efficiencies are 0.4(1) 
and 5.2(4) for 87 Rb and 133 Cs respectively, as we sacrifice 
87 Rb to sympathetically cool 133 Cs (□ and ■ in Fig. [2] 
(i)). As the dipole trap is loaded the RF frequency is re- 
duced to optimise the transfer of 133 Cs whilst at the same 
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FIG. 2: (Color online) Trajectory to 133 Cs BEC via sympa- 
thetic cooling with 87 Rb (squares). Open (closed) symbols 
show 87 Rb ( 133 Cs) data. Circles show single species data for 
comparison. The trajectory gradient indicates the cooling ef- 
ficiency. The different regions show: (i) RF evaporation in the 
magnetic quadrupole trap, (ii) Dipole trap loading and inter- 
nal state transfer, (iii) Evaporation by reducing the dipole 
trap beam powers, (iv) Evaporation by trap tilting. Inset: 
3D rendering of the density distribution as the cloud is con- 
densed. Right: absorption images (350 x 350 /xm) and cross 
sections of the cloud near the BEC phase transition. 



time removing all of the 87 Rb. The resulting loading is 
highly efficient with ~ 50 % of the 133 Cs atoms trans- 
ferred into the dipole trap at a temperature of 10(1) fiK. 
The phase-space density is 6.6(2) x 10 -4 , representing 
more than a factor of 400 increase compared to the final 
value in the quadrupole trap. This increase results from 
an effect analogous to the use of a dimple potential [HI . 
Adiabatic expansion of the trap volume cools the cloud 
in the quadrupole potential and elastic collisions transfer 
some atoms into the tighter harmonic optical potential 
[35| ] . The remaining atoms are lost from the quadrupole 
trap reservoir, removing energy in the process. We be- 
lieve the enhanced loading of 133 Cs in the presence of 
87 Rb compared to the 133 Cs alone case (• in Fig. [2j) is 
again consistent with the existence of a large interspecies 
scattering length [33|, [34]. 

To proceed we transfer the 133 Cs atoms to the |3, +3) 
state with ~ 95 % efficiency using RF adiabatic rapid 
passage as a bias field is turned on to ~ 22 G in 18 ms 
[30|. This transfer eliminates inelastic two-body losses 
which have plagued previous attempts to cool 133 Cs to 
degeneracy in magnetic traps 136]. In addition, at this 
bias field the ratio of elastic to inelastic three-body colli- 
sions is favorable for evaporative cooling [37| . Operating 
in the vicinity of the broad zero-crossing in the scatter- 
ing length near 17 G gives precise control of the elastic 
collision rate. We observe efficient evaporation for fields 
between 21 and 24 G, and typically operate at 22.4 G 
where ac s = 280 ao- The majority of the evaporation is 
performed by reducing the beam powers to 48 mW over 
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FIG. 3: (a) Interspecies inelastic losses in the dipole trap. 
To determine trap lifetimes the majority species ( 87 Rb) data 
are fit with a double exponential function while the minority 
species ( 133 Cs) data shown in the inset are fit with a single 
exponential function, (b) Trajectory to a dual-species BEC 
of 87 Rb and 133 Cs. The trajectory is divided into four sec- 
tions: (i) RF evaporation in the magnetic quadrupole trap, 
(ii) Dipole trap loading and internal state transfer, (iii) Evap- 
oration by reducing the dipole trap beam powers, (iv) Evap- 
oration by trap tilting. Open (closed) symbols show data for 
87 Rb ( 133 Cs). 



9.5 s. The final evaporation to degeneracy is implemented 
by increasing the magnetic field gradient to 35.4(1) G/cm 
over 2.0 s. This tilts the trapping potential (Fig.[T] (d)), 
thereby lowering the trap depth whilst leaving the trap 
frequencies largely unchanged (29j . With this approach 
we achieve an evaporation efficiency of 2.1(1) for 133 Cs 
alone in the dipole trap and produce pure single species 
condensates of up to 6.2(1) x 10 4 atoms. Fig. [2] sum- 
marizes the trajectory to 133 Cs BEC. We are not able to 
produce condensates by loading only 133 Cs into the mag- 
netic trap due to the poor transfer into the dipole trap 
(• in Fig. [2j). In contrast, loading only 87 Rb initially we 
are able to form condensates of up to 1 x 10 6 atoms (30| . 

To create a degenerate mixture of 87 Rb and 133 Cs we 
reduce the number of 133 Cs atoms loaded into the mag- 
netic trap to 2.7(1) x 10 6 . This reduces the heat load on 
the 87 Rb during the sympathetic cooling and improves 
the cooling efficiencies in the magnetic trap to 1.7(1) for 
87 Rb and 11.3(4) for 133 Cs. Loading the dipole trap re- 
sults in a mixture of 2.8(2) x 10 6 87 Rb and 5.1(3) x 10 5 
133 Cs atoms at a temperature of 9.6(1) /iK following the 
spin flip to the 87 Rb |1, +1) and 133 Cs |3, +3) states . The 
initial peak densities for both species are > 10 13 cm -3 
and we observe very strong interspecies inelastic losses, 
see Fig. [3fa). The fast loss of 133 Cs atoms from the 
trap shown in the inset of Fig. [3]^a) is purely due to 
inelastic three-body collisions. The initial decrease of 
the number of 87 Rb atoms includes a significant con- 
tribution from plain evaporation from the trap. Three 
body collisions are most prevalent at the trap centre and 
cause the coldest atoms to be lost predominantly. This 
'anti-evaporation' leads to both atom loss and heating 
and is a big obstacle when trying to produce a dual- 



species BEC. The heating leads to further loss from the 
trap through evaporation and complicates the determi- 
nation of the three-body loss rate. A full analysis of 
the evolving densities using two coupled differential rate 
equations [38| is therefore not performed. Instead we 
estimate an upper limit for the interspecies loss rate co- 
efficient of ~ 10 -25 — 10 -26 cm 6 s _1 based simply upon 
the observed lifetime of the minority species and the av- 
erage initial densities. The three-body loss rate coeffi- 
cient scales as L3 = 3C— a 4 , where C is a dimension- 
less factor between and ~70, m is the reduced mass 
and a is the scattering length [39|. The current best 
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34j of the interspecies scattering length, 
650 ao, would lead to a maximum three-body 
loss rate of L3 ~ 3 x 10 -25 cm 6 /s in reasonable agreement 
with our simple estimates. 

To combat the strong interspecies inelastic losses the 
beam powers are reduced to 120 mW in just 1.0 s, cor- 
responding to a final trap depth of 2 /iK for 87 Rb. Tak- 
ing into account the adiabatic expansion of the cloud 
resulting from the commensurate relaxation of the trap 
frequencies, this corresponds to an effective reduction in 
the trap depth by a factor of 7. The result is a factor 
of ~ 3 reduction in the density and consequently an or- 
der of magnitude increase in the lifetime of 133 Cs against 
interspecies three-body collisions. Overall the conditions 
for evaporative cooling are improved by this fast reduc- 
tion in the dipole trap power, despite the associated de- 
crease in the elastic collision rate. Further evaporation is 
then performed by reducing the beam powers linearly to 
48 mW in 2.5 s. Finally by again tilting the trap using the 
applied magnetic field gradient dual-species condensates 
are produced in the same trapping potential containing 
up to ~ 2.0 x 10 4 atoms of each species. Similar results 
have recently been obtained by confining the 87 Rb and 
133 Cs atoms in spatially separated optical traps to sup- 
press the interspecies inelastic losses [34|. The trajectory 
to a dual-species BEC is presented in Fig. [3^b). The 
efficiencies of the evaporative and sympathetic cooling in 
the dipole trap are 1.8(2) and 2.8(5) for 87 Rb and 133 Cs 
respectively. The data show that despite strong inter- 
species inelastic losses the sympathetic cooling of 133 Cs 
via 87 Rb in the dipole trap is still more efficient than the 
direct evaporation of 133 Cs alone (Fig. [2]). 

A dramatic spatial separation of the two condensates 
in the trap is observed, revealing the mixture to be im- 
miscible at 22.4 G (Fig. |4]). The resulting density profile 
for each species is a stark contrast to the symmetry of 
the potential and the density profiles observed for each 
single species condensate. For an immiscible quantum 
degenerate mixture, the relative strength of the atomic 
interactions A = ^ RbC , s . > 1 , where the interaction cou- 
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&Rb = 100 ao and ac s = 280 ao, so that the observation 
of immiscibility implies aRbCs > 165 ao; again consistent 
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FIG. 4: (Color online) Immiscibility in a quantum degener- 
ate mixture of 87 Rb and 133 Cs. Three distinct structures are 
observed (A, □ and •) correlated with the number of atoms 
in each condensate. Absorption images and cross-sections of 
87 Rb (red) and 133 Cs (blue) highlight these structures. The 
optical depth (OD) scales of the cross-sections are 0-1.0 (a) 
and 0-2.6 (□ and •). Each image is 350 x 100 fim. 
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with the current best estimate, aRbCs ~ 650 ao 
The bimodal distributions in Fig. |4]show that the immis- 
cible behavior is exclusive to the condensed atoms. The 
immiscible condensates always form one of three struc- 
tures, either one of two possible symmetric cases (a and 
• in Fig. 3]) or an asymmetric case (□ in Fig. 3]). Qual- 
itatively, these observations are strikingly similar to the 
symmetric and asymmetric structures predicted for an 
immiscible binary condensate based upon numerical sim- 
ulations of two coupled Gross-Pitaevskii equations [41 1. 
Experimentally, the structures are correlated with the 
number of atoms in each condensate. Altering the num- 
ber of 133 Cs atoms loaded into the trap, allows navigation 
between the three structures. The strong interspecies 
losses make the formation of a dual-species condensate 
with > 2 x 10 4 atoms of each species a challenge and 
lead to the asymptotic nature of the data in Fig. 2J 

In summary, we have demonstrated the production of a 
dual-species Bose-Einstein condensate of 8T Rb and 133 Cs 
in the same trapping potential. Our method uses simple 
established techniques and exploits the large interspecies 
scattering length to efficiently sympathetically cool 133 Cs 
in both the magnetic and optical trapping phases of the 
experiment. We observe a striking phase separation of 
the two condensates revealing the mixture to be im- 
miscible due to the strong repulsive interspecies inter- 
actions. We believe that the immiscibility aids the for- 
mation of the dual-species condensate by suppressing the 
interspecies inelastic losses as the mixture condenses. In 
the future, we plan to explore the miscibility by Feshbach 
tuning both ac s and aRbCs m order to test this assertion. 



Crucially, thermal clouds very close to degeneracy remain 
miscible, enabling the efficient magneto- association [21] 
of RbCs molecules in future experiments, as the first step 
towards the creation of a dipolar molecular quantum gas. 
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